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This is an analytical study of pinning and spontaneous vortex phase is a system consisting of a 
superconducting thin film pierced by a long ferromagnetic columnar defect of finite radius R. The 
magnetic fields, screening currents, energy and pinning forces for this system are calculated. The 
interaction between the magnetic field of vortices and the magnetization outside the plane of the film 
and its close proximity enhances vortex pinning significantly. Spontaneous vortex phase appears 
when the magnetization of the columnar defect is increased above a critical value. Transitions 
between phases characterized by different number of fiux quanta are also studied. These results are 
generalized to the case when the superconductor is pierced by an array of columnar defects. 

PACS numbers: 74.25.Ha, 74.25.Qt, 74.78.-w, 74.78.Na 



I. INTRODUCTION 

Optimizing pinning in superconductors (SC) is a prob- 
lem that is interesting both experimentally and theoret- 
ically since pinning is the main viable mechanism for 
superconductivity in the presence of external magnetic 
field. Many tools and mechanisms for pinning has already 
been studied, in particular the use of crystal defects such 
as holes, non-magnetic impurities and both linear and 
screw dislocations. Pinning using structural defects suf- 
fers from many drawbacks, the most important is that 
the pins are randomly distributed which results in a low 
critical current. In recent years, it was claimed that pin- 
ning could be optimized if we employ ferromagnetic (FM) 
nano textures to pin superconducting vortices (SV)»i^ 

The main obstacles of using ferromagnetic textures 
to pin superconducting vortices are the proximity ef- 
fects which destroy superconductivity. Lyuksyutov and 
Pokrovsky noticed that proximity effects could be sup- 
pressed if a thin layer of insulator oxide is sandwiched 
between the FM and the SC.^ Due to this separation be- 
tween the FM and SC, the interaction between them is 
mediated via their magnetic fields. If the magnetization 
of the FM structure exceeds a threshold value Mc, the 
interaction between the FM and SC makes the sponta- 
neous creation of superconducting vortices energy favor- 
able. The interaction between the FM and SC increases 
as the temperature reaches the superconducting transi- 
tion temperature Tg from below. 

Marmorkos et. al& considered a system consisting of a 
thin ferromagnetic dot embedded in a superconducting 
thin film. Both the superconductor and the ferromagnet 
are assumed to lie in the sy-plane and the dot is magne- 
tized along z-axis. By numerically solving the nonlinear 
Ginzburg-Landau equation, they were able to show that 
a vortex appears in the superconductor when the mag- 
netization of the dot exceeds a critical value. They also 
showed that increasing the dot's magnetization leads to 
a giant vortex state with multiple flux quanta. 

In Ref0, Erdin et.al. used London's theory of super- 
conductors to solve the general problem of the interac- 



FIG. 1: A superconducting thin film pierced by a ferromag- 
netic nano rod of radius R, length 2L and magnetization M. 



tion between vortices in superconducting thin film with 
a generic two-dimensional ferromagnetic structure. They 
applied their results to the cases when the FM structure 
is a circular dot whose magnetization is either parallel or 
perpendicular to the plane of the superconductor. They 
calculated the threshold value for the dot's magnetization 
at which vortices are spontaneously created in the super- 
conductor. They showed that by increasing the magne- 
tization of the dot a series of phase transitions between 
phases with different number of vortices take place. 

All these studies focused on cases in which the ferro- 
magnet is either a point-like dipole or an infinitely thin 
two-dimensional texture whose plane is parallel to the 
plane of the superconductor. The problem of interaction 
between vortices in superconducting thin film with a fer- 
romagnetic columnar defect (FCD) has not been studied 
yet. This case is both theoretically and experimentally 
interesting for the following reason. The magnetic field of 
the vortex in a superconducting thin film is not confined 
to the plane of the film but it also exist outside of it. 
If the ferromagnet extends in space outside the film and 
its close proximity then the vortex will be bound more 
strongly to the ferromagnet due to the interaction of its 
magnetic field with the magnetization of the ferromag- 
net. This enhances vortex pinning and consequently the 
critical current of the superconductor. It is also impor- 
tant to note that ferromagnetic columnar defects makes 
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the superconductor has a multiply connected topology 
which is another source of pinning. Therefore, it would 
be interesting both theoretically and experimentally to 
study the statical and dynamical properties of a super- 
conducting thin film pierced by ferromagnetic columnar 
defects. 

In this article, I propose to study the properties of 
spontaneous vortex phase and pinning in a superconduct- 
ing thin film pierced by ferromagnetic columnar defects. 
This article is organized such as in the first section, the 
magnetic potential and field distribution produced by a 
single FCD penetrating a thin superconducting film are 
calculated. In sections two, I calculate the total energy 
of a system of a superconducting vortex coupled to the 
magnetic defect. I generalize these results to the case 
of a superconducting thin film pierced by a square ar- 
ray of magnetic columnar defects in section four then I 
summarize the results of this work in section five. 



II. AN SC FILM PIERCED BY A 
FERROMAGNETIC COLUMNAR DEFECT 

I consider a superconducting thin film in the xy-plane 
pierced by a finite ferromagnetic nano rod whose radius 
is R and length is L as shown in Fig.Q] I also assume that 
the FCD is uniformly magnetized along its symmetry axis 
then its magnetization distribution is written as 



The total current J = J** -I- J™ is the sum of the super- 
current J" and the magnetization current J™ = cV x M. 
The supercurrent can be written in terms of the gradient 
of the superconducting order parameter x and the total 
magnetic vector potential A = A^ + A™ for the FCD-SC 

system J" = §^ (vx - f^A^ . Here is the Super- 
conducting electron density, h is Planck's constant while 
TOg and e are the electron's mass and charge respectively. 
The superposition principle allows us to solve (O for A™ 
and A* independently. If the Coulomb gauge V • A™ = 
is imposed then A™ satisfies the following equation 



V^A™ = -iA™(5(0) + 47rV x M 
A 



(3) 



where [A = Xj^/dg] with [Al = y^mgC^ / iirnse'^] is Lon- 
don's penetration depth. Here rUe and e are the electron's 
mass and its charge while is the superconducting elec- 
trons density and dg is the thickness of the SC film. The 
general solution of Q is of the form 



A™ = J A'^e-<'i-P+''''^d^qdk, 



(4) 



where and q are the components of the wave vector 
along the 2;-axis and in the xy-p\ane. Here is the 
vector potential in momentum space and is given by 



MeiR-p)ei--\z\ 



(1) 



where z is the unit vector along the 2:-axis. The magnetic 
field produced by the FCD penetrates the superconduc- 
tor and changes the distribution of the screening current. 
In the London approximation the FCD-SC system is de- 
scribed by London-Maxwell equation 
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V X V X A = — J 
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where Jn{x) is the n-th rank Bessel's function. The 
magnetic field produced by the magnetic nano-rod in 
the presence of the SC film can be calculated using 
B = V X A™. The components of the rod's magnetic 
field are 



SnMR 
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where W{q, z, L) is 

W{q,z,L) — sign{L — 2z) 

-\- sign{L -\- 2z) 



^W{q,z,L) 



^Wiq,z,L) 



sinh(^)e-9(l^l+*) 
l + 2q 

sinh(^)e-^(l"l+^) 
l + 2q 



dq 
dq 



, ,q(L - 2z), ,^ „ X . -, ,qiL - 2z) ' 
1 - cosh(-i^^ ^) + sign{L - 2z) smh{-^ '-) 



,,q{L + 2z). . , .q\L + 2z\ 
1 - cosh(-i^^ — -) + s.\n\i{- — 



(6) 
(7) 



(8) 



However, we are interested in the value of the field at the 
plane of the superconductor. The z component of the 



magnetic field of the FCD evaluated at the SC plane is 

BTip)^'-^ r^:^lM^(l-e-^)d, (9) 
A Jo 1 + 2g 
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Similarly the solution for A* is found. A general argu- 
ment made ini shows that the term proportional to Vx 
ascribes for vortices. So, in the presence of a vortex with 
vorticity v and center at po = 0, the solution for A^ in 
the Coulomb gauge is 



■/i(g|p|)e"''l^l 
l + 2Ag 



dq (10) 



where p is the unit vector along p. The z-component of 
the vortex magnetic field at the SC film surfac eAi.. ji 5 i is 



27rA2 



Yp 



(11) 



where Hq{x) and iVo(a;) are the zero order Struve and 
Neumann functions respectivelj>i&. 



III. THE PINNING POTENTIAL AND ENERGY 
OF FCD-SV SYSTEM 

The interaction between a superconducting vortex and 
a non magnetic columnar defect was first considered 
by Mkrtchyan and SchmidUi and later in the work of 
Buzdin et.alm^^. In these studies, it was shown that the 
pinning potential Up created by a non magnetic columnar 
defect of radius R > is 



Upip) 



-^oln(^), 



eo In 



1 _ ( 



p<R 
R < p < X 



where ^ is the SC coherence length and eo — is the 

energy scale of the vortex self interaction. 

If the columnar defect is ferromagnetic, then an extra 
contribution to the pinning would appear due to the in- 
teraction between the superconductor and ferromagnet. 
In the presence of vortices in the superconductor, the 
total energy of the system is made up of five different 
contributions and can be written as 



U = Usv + Uvv +Up + Umv + Ur, 



(12) 



where Usv is the energy of N non-interacting vortices, 
U^^ is the vortex-vortex interaction, Umv is the interac- 
tion energy between the FM and the SC, and Umm is the 
FM dot self interaction. In^, it was shown that the total 
energy of the system can be written as: 



U = 



8m 



(Vx • A) - -M • B 

Arrive 1 



where c is the speed of light. The vectorial quantities 
A, and B are the total vector potential and magnetic 
field due to the vortices and the ferromagnetic columnar 
defect. 



In Refl23l Buzdin considered the problem of vortex 
craetion in a superconductor penetrated by non-magnetic 
columnar defects and placed in external magnetic field. 
He found that at high magnetic field giant vortices are 
more energy favorable when the columnar defect is rather 
thick. The case of rather thick defects i?/A 1 is exper- 
imentally more interesting and easier to realize. Mar- 
morkos et.aL? showed the possibility of formation of a 
giant vortex with multiple flux quanta around the mag- 
netic dot instead of singly quantized vortices. Therefore, 
I will assume that ^ <C i? < A which make the result of 
Reffeil regarding the formation of giant vortices applica- 
ble to our problem. 

The phase gradient of the SC order parameter in the 
presence of a giant superconducting vortex with vortic- 
ity 1/ is Vx — t^ ^^pl^p^p^ ; where po is the location of 
the vortex. The total energy for a system consisting of 
a superconducting vortex and a ferromagnetic columnar 
defect is 



Uipo) 



v^eo In(-) + i/2t/p(po) 



(14) 



^ JO 



Jo(<zpo)Ji(<zi?)(l-e-V) 



dq 



g(l + 2Aq) 

where — A/^qA and Umm is the self interaction of 
the FCD is ignored since it does not affect the super- 
conducting state. It is clear that if the magnetization of 
the columnar defect is large enough so that the interac- 
tion term becomes larger than other terms in ((Tljl then 
the vortex can appear spontaneously in the supercon- 
ductor. Erdin et.ali showed that spontaneous creation 
of vortices takes place most easily in the proximity of the 
superconducting transition temperature T,. 

The interaction term in Ea. (|14|l depends on L. In the 
limit of L/R ^ 1, the interaction term in Ea. H14|l reduces 
to the result obtained for infinitely thin magnetic doti 
In the remainder of this work, I will assume L oo. 
Numerical minimization of the total energy U sets pa = Q, 
hence the vortex center must be on the axis of the FCD. 

The energy of a system of a singly quantized vortex 
V = \ coupled to a columnar defect is shown in Fig. |2] 
The solid line represents the energy of the system when 
the columnar defect is non-magnetic, while the dashed 
line shows the energy of the system when the defect is 
ferromagnetic. Note that the pinning force, — VJ7(p), is 
not zero for p < i? in contrast with the case of non- 
magnetic columnar defect. 

For each value of M there is a corresponding critical 
value of v which I call it here Vc- The minimization of 
C/(po = 0) with respect to u yields 



167rmo R r°° JijqR) 

ln(A) A X q{i + ^M) 



dq 



(15) 



where mo = M/Mq with Mq = (^^o/ttA^. The critical 
value of v is the closest integer to the value of ly given by 
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FIG. 2: The total energy of the CD-SV system for A = 
lOOOnm, ^ — lOnm and the radius R — 900nm.. The sohd 
hne is for the case when the CD is non-magnetic while the 
dashed one is for ferromagnetic CD. 
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FIG. 4: The curve represents the threshold value of the mag- 
netization of the FCD in units of Mo = (^o/ttA^. All points 
above the curve represent a phase with vortices while the re- 
gion under the curve is vortexless. 



Ea. (|15|l . The vorticity of the spontaneously created giant 
vortex is plotted in Fig. 13 as a function of toq and R/X. 
Note that at when R/X ^ 1 then very large values of 
Too are required for the spontaneous creation of a vortex. 
Note that for any fixed value of i^c, Ea. H15() gives us the 




FIG. 
R/X. 
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3: (Color Online) The plot of j/ as a function of mo and 



while its energy of interaction with the columnar defect 
is not zero and negative; therefore, giant vortices with 
large value of u are energy favorable. 

The curve in Fig. 21 is the curve for the first phase 
transition (No vortex Phase — > One vortex phase), it 
should be emphasized here that there is a discrete set 
of curves above this curve that refer to a series of phase 
transitions between phases with different number of flux 
quanta. The transition from a vortexless phase to a phase 
with a single vortex for a system in which R — 50 nm 
and A — 100 nm occurs at the approximate value of 
Mc = 66.85 G. 



IV. AN SC FILM PIERCED BY A SQUARE 
ARRAY OF FCDS 



Now, let us consider a superconducting thin film 
pierced by a square array of N x N ferromagnetic colum- 
nar defects. Let each FCD be infinitely long {L oo) 
with a radius ^ <C i? < A and magnetization AI > 
Mcivc)- The magnetization of the array can be written 
in Fourier space as follows 



critical value of the magnetization Md^c) at which a vor- 
tex with vorticity Vc appear as a function of the radius R. 
The dependance of Mc{vc) on R for Vc — 1 \s depicted in 
Fig. 21 The region under the curve in Fig. ^represents a 
vortexless phase while the one above the curve represents 
phases with vortices. 

Note that when the radius of the columnar defect is 
very close to A giant vortices with large flux quanta are 
expected even at not so large values of toq. This is be- 
cause when i? « A, the vortex self energy goes to zero 



M 



K = AttHIR^^^5{K)z y e'^ P^ (16) 
1 ~ 



where pi is the location of the «-th FM defect and the 
sum runs over all the lattice sites. The location of the 
FCD is defined by two integers I and s as pi = lax + say 
where a is the lattice spacing of the FCD array. Solving 
Eq.(|3Jl for the magnetization H16|) yields 
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\&Ti'^iMRJi{qR) 
X (q X z) 
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q{l + 2\q) 



Pj 



(17) 



Let Q be a vector of the reciprocal lattice defined as 
Q = TiQi + r7iQ2 where n and m are two integers with 
Qi = and Q2 = "^V- The z component of the 
array's magnetic field at z = is 

B^{x,y,z^ 0) = 16^2^oMAi?^ A^lge-Q '' (18) 



where no = 1/a^ is the density of colmimar defects in 
the array. The summation over the defect positions is 
calculated with the help of the identity 



(2^)2 



Q 



(19) 



The magnetic fiux which penetrate the superconductor 
must remain constant; therefore, the number of vortices 
and the number of antivortices that are spontaneously 
created in the superconductor must be equal. In the ab- 
sence of external magnetic field, the equilibrium configu- 
ration is the one with vortices pinned at the axes of the 



FCDs and antivortices are located at the centers of the 
unit cells of the FCD lattice. The square lattice of an- 
tivortices is identical to that of the vortices with a shift 
by a/2 — {a/2, a/2). Note that in the case of a single 
columnar defect, the antivortex is repelled (pushed far 
away) from the FCD; therefore, it was not taken into ac- 
count in the calculations. However, in the case when the 
SC is penetrated by an array of FCD in the absence of 
external magnetic field, one must take antivortices into 
account on equal footing with vortices. In equilibrium, 
the force acting on any flux line (FL) in the supercon- 
ductor is zero. Pinning forces appear as soon as this 
equilibrium is disturbed e.g. by passing an electric cur- 
rent to the superconducting film. The remaining part of 
this section will focus on the problem of calculating these 
forces. The vortex and antivortex lattices are both regu- 
lar and periodic. Therefore, if a small current is applied 
to the superconductor then all vortices will move together 
in one direction and all antivortices will move together 
in the opposite direction. Let us assume that the vortex 
lattice is displaced with respect to its equilibrium posi- 
tion by a small amount Ap = (cc, y) then the lattice of 
antivortices will be displaced by an amount — Ap. Con- 
sidering the case when M is just above M(,{vc = 1). The 
fiux line self energy and the FL-FL interaction energy 
together can be written using Ea. H13|l as 



Anl(j)l 



?Q(1 



2t: ^Q(1 + 2AQ) 



l-cos(Q.(--2Ap)) 



(20) 



where A is the area of the SC film. The energy of interaction between the flux lines and the ferromagnetic array is 

MQR) 



E^, = -27r AnlMRXS2-—^ 
q Q (1 



Q2(1 + 2AQ) 



cos(Q • Ap) — cos(Q • (— — Ap) 



(21) 



r 



To simplify the calculations, I use 



(22) 



where ^0 = is the area of the array's unit cell. There- 
fore, I find the energy per flux line is now dependent on 
the displacement Ap as follows 



p - A. 

° " 87r2 



[1~ Jo((z|f-2Ap|)] 



dq 



-M(j)oR 



l + 2Xq 

JijqR) [Jo(g|Ap|)- Jo(g||-Ap|)] 
q{l + 2\q) 



(23) 

dq 



Note that the number of FCD is N"^ while the number of 
flux lines is 2N'^ . The pinning force acting on any vortex 



or antivortex can now be calculated using fp — —VU{Ap) 
where ?7(Ap) = Eo — £e with is a constant equals to 
the energy of the system in equilibrium. Note that I did 
not include the term due to the Mkrtchyan and Schmidt 
in the result I obtained in Ea. H24|) because it is rather 
straightforward to calculate its effect. 



V. CONCLUSION 

In conclusion, I studied the interaction between super- 
conducting vortices in a thin film pierced by a ferromag- 
netic columnar defect. I calculated the magnetic field of 
the FCD in the presence of the SC and the distribution 
of screening currents in the superconductor. If the mag- 
netization of the FCD exceeds a critical value, then the 
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FCD interaction with the vortex will overcome the vortex 
self energy leading to the spontaneous creation of vortices 
in the superconductor. I showed that vortex pinning is 
strongly enhanced due to the contribution from the in- 
teraction between the magnetic field of vortices and the 
magnetization of the FCD outside the plane of the SC 
and its close proximity. This extra contribution to vor- 
tex pinning is a major difference between this study and 
other studies in the literature which focused on dipole- 
like or two-dimensional ferromagnetic structures. These 
results were generalized to include the case when the SC 



film is pierced by an array of FCD. A detailed analysis of 
the dynamical properties of this system will be reported 
elsewhere. 
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